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ABSTRACT 


The  problem  of  determining  the  optimum  allocation  of  aircraft  to 
an  alrstrike  against  a  transportation  network  is  investigated.  The 
damage  function  is  assumed  to  be  exponential.  A  solution  procedure  is 
developed  utilizing  dynamic  programming  and  integer  solutions  are  found. 
The  number  of  aircraft  to  be  assigned  to  the  airstrike  is  considered  a 
decision  variable.  A  sensitivity  analysis  is  run  to  determine  the 
optimum  value  for  this  variable. 
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I.  INTRODUCTION 


A.  OBJECTIVE 

The  purpose  of  this  paper  is  to  present  a  procedure  for  determining 
the  optimal  allocation  of  aircraft  to  a  single  airstrike  against  a 
transportation  network.  This  allocation  problem  is  solved  by  dynamic 
programming  and  a  fortran-coded  version  of  the  program  is  included  in 
the  paper. 

B.  GENERAL 

Sustained  ground  operations  require  a  military  force  to  have  some 
means  of  resupply.  This  resupply  capability  is  partially  dependent 
upon  a  land  transportation  system.  The  level  of  resupply  effort  required 
depends  upon  what  type  of  forces  are  being  supported.  Guerrilla  forces 
enjoying  local  support  require  less  resupply  capability  in  terms  of 
pounds  per  man  per  day  than  would  a  conventional  army,  but  a  greater 
percentage  of  this  capability  depends  upon  land  transportation  networks. 

Any  reduction  in  the  resupply  capability  of  a  military  force  will 
reduce  its  combat  effectiveness.  Tactical  air  interdiction  has  been 
used  extensively  by  the  Armed  Forces  of  the  United  States  against  its 
opponents  in  Southeast  Asia  to  accomplish  this  reduction. 

There  are  at  least  three  alternative  means  of  using  tactical  air  to 
reduce  the  resupply  capability  of  an  enemy.  Aircraft  may  be  assigned 
to  attack  sources  of  supply  to  destroy  war  material  before  it  enters 
the  transportation  system  and/or  to  disrupt  its  production;  aircraft 
can  destroy  war  material  as  it  moves  in  the  transportation  system;  and 
finally  aircraft  can  attempt  to  reduce  the  resupply  capacity  of  the 
transportation  system  itself  by  destroying  bridges,  roads,  railroads, 
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et  cetera.  Conventional  wisdom  argues  that  the  first  course  of  action 
Is  the  most  effective  form  of  interdiction.  Unfortunately  for  military 
planners,  political  considerations  may  rule  out  this  alternative.  This 
paper  will  focus  on  the  last  of  these  options,  the  reduction  in  capacity 
of  the  transportation  system  itself. 

C.  BACKGROUND 

Considerable  effort  has  been  devoted  to  the  interdiction  problem. 

In  particular  two  recent  papers  provided  the  background  for  the  approach 
to  the  problem  developed  in  this  paper.  McMasters  and  Kustin  [1] 
developed  an  algorithm  that  determines  which  arcs  of  a  transportation 
network  should  be  attacked  and  at  what  level  of  effort  given  a  limited 
availability  of  resources.  In  this  formulation  of  the  problem  the 
relationship  between  arc  capacity  and  resource  allocation  (damage 
function)  was  assumed  to  be  linear.  The  algorithm  presented  is  based 
upon  the  max-flow  min-cut  theorem  of  Ford  and  Fulkerson  [2]  and  the 
relationship  between  a  primal  network  and  its  topological  dual. 

Nugent  [3]  investigated  the  same  problem  under  the  assumption  of  an 
exponential  damage  function  which  exhibits  diminishing  marginal  returns. 
An  algorithm  was  developed  that  finds  a  non-integer  solution  to  the 
problem. 

In  this  paper  the  transportation  system  will  have  the  same  network 
formulation  as  in  Refs.  1  and  3.  The  problem  will  be  formulated  differ¬ 
ently  and  dynamic  programming  will  be  used  to  provide  integer  solutions. 

D.  INTERDICTION  PROBLEM 

It  will  be  assumed  that,  given  unlimited  aircraft  availability,  the 
assignment  of  aircraft  to  an  airstrike  would  reach  a  point  beyond  which 
it  would  become  uneconomical  to  assign  further  aircraft.  In  a  problem 
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with  constraints  on  aircraft  availability  this  point  might  or  might  not 
occur  before  all  available  aircraft  were  assigned.  For  this  reason, 
the  objective  of  an  operations  officer  planning  an  airstrike  against  a 
transportation  network  is  not  merely  to  minimize  network  capacity 
subject  to  aircraft  availability,  but  to  minimize  the  capacity  subject 
to  aircraft  availability  and  the  additional  consideration  that  the  cost 
of  any  Incremental  assignment  of  aircraft  to  the  strike  is  exceeded  by 
the  benefit  resulting  from  that  assignment. 

To  accomplish  the  objective  the  strike  planner  must  have  information 
on  the  availability  and  cost  of  assignment  of  aircraft.  Detailed  infor¬ 
mation  must  be  available  concerning  the  transportation  network  including 
the  upper  and  lower  bounds  on  the  capacity  of  each  arc  and  its  vulnera¬ 
bility  to  attack.  The  planner  must  also  know  the  benefit  to  attribute 
to  a  reduction  in  resupply  capability.  With  this  information  and  using 
the  procedure  that  will  be  outlined  the  planner  can  determine:  how 
many  aircraft  to  assigne  to  the  airstrike;  which  arcs  in  the  network 
should  be  attacked;  how  many  aircraft  to  as  ign  to  arcs  that  will  be 
attacked;  and  the  capacity  of  the  network  after  the  airstrike. 


II.  THE  MODEL 


A.  NETWORK  DESCRIPTION 

The  transportation  system  under  consideration  is  represented  by  a 
planar  connected  graph  of  nodes  and  undirected  capacitated  arcs.  Arcs 
represent  road  segments  and  nodes  represent  either  a  road  intersection 
or  any  other  point  where  it  is  necessary  to  distinguish  between  road 
characteristics  on  either  side  of  the  node.  Three  constants  are  asso¬ 
ciated  with  each  arc  representing  the  upper  and  lower  bounds  on  arc 
capacity  and  the  arc's  vulnerability  parameter. 

It  is  assumed  that  the  network  has  one  source  node  from  which  flow 
originates  and  one  sink  node  at  which  flow  terminates.  If  the  trans¬ 
portation  system  being  modeled  has  more  than  one  originating  point  or 
terminating  point  this  may  be  handled  by  creating  a  super-source  and/or 
sink  with  artificial  arcs  connecting  these  super-nodes  to  sources  and 
sinks  as  needed.  These  artificial  arcs  may  not  be  attacked  and  their 
capacities  are  unbounded.  The  arc  between  nodes  i  and  j  is  represented 
by  (1,j).  Nodes  are  numbered  from  1  to  n  with  1  corresponding  to  the 
source  and  n  the  sink.  With  the  exception  of  the  source  and  the  sink, 
flow  conservation  is  assumed  to  hold.  That  is,  flow  out  of  node  i 
equals  flow  Into  node  1. 

The  flow  in  arc  (1,j)  Is  designated  as  if  It  is  from  node  1  to  j 
and  Xjj  If  it  is  from  node  j  to  i.  This  avoids  the  necessity  of  defining 
negative  flows.  Flow  Is  assumed  to  be  from  the  source  to  the  sink 
although  It  may  be  In  either  direction  In  the  Intermediate  arcs.  The 
model  as  formulated  considers  only  flows  of  a  single  commodity,  tons 
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of  resupply  per  day.  ard  the  value  of  one  unit  of  flow  is  assumed  to 
be  the  same  for  aP  arcs. 

Capacities  on  arcs  represent  bounds  cn  flow  In  either  direction. 

The  capacity  on  arc  (1,j)  is  given  by  riy  and  is  assumed  to  be  the  same 
In  both  directions.  The  flow  in  arc  (i,j)  is  restricted  by 


-"Hj  • 

The  upper  and  lower  hour..*  on  the  capacity  of  arc  (1,j)  are  repre¬ 
sented  by  Uy  and  ly  where 
* 

Ol'ij  • 


The  vulnerable  portion  of  an  ire's  capacity  is  designated  Wy  with 


WU  a  utj  - 

The  amount  of  resource  allocated  to  interdict  arc  u,j)  is  denoted  by 
ky.  The  relationship  between  the  .apecity  o;  arc  1 1 . j )  and  the  level 
of  resource  assigned  to  its  nrerctetion  defined  as  "he  danage 
function  of  arc  (i,j)  and  is  given  by 


■  >lj  *  \j  s"?(-ouk,.)  . 

In  the  above  damage  funcJcn  tru  parameter  by  ij  a  measure  of  the 

vulnerability  of  arc  (i,j).  La*  er  voices  of  by  result  in  greater 

reductions  in  capacity  for  fi  ed  'u'oes  of  1^,  *y  and  ky  and  hence 

Imply  greater  vclntraoil  ity.  It  i»y  •  r'  ’.her  my>  ^ )  «  uy  for 

all  possible  values  of  kj,  and  no  cencity  •  -’duction  i  possible  With 
this  damage  function  if  no  eircra*t  are  assigned  ti  (1,i),  4ts  capacity 
will  be  uy  and  in  the  Hint  as  the  numre-  of  aircraft  assNned  to  (i,j) 
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becomes  Infinite  tne  capacity  approaches  This  lower  Lojnd  Kill 

be  referred  to  as  arc  capacity  after  unlimited  interdiction. 

B.  DETERMINATION  OF  NETWORK  CAPACITY 

The  opposition  is  assumed  to  have  the  oeans  to  oetermire  how  to 
maximize  the  flow  in  the  transportation  network.  Let  tne  capacity  of 
the  network  be  defined  as  this  maximal  flow,  i he  detennnatic.  of 
maximum  flow  Is  the  well-kncwn  maximal  flow  problem  ano  -wy  be  found 
using  the  max-flow  lac.*'. ing  algorithm  base*  jpon  the  max- flow  m in-cut 
theorem  of  Tcra  ana  cu  Ike*  son  [Z].  Ford  and  Fulkerson's  theoro.i  states 
that  the  maximu.il  flow  possible  in  e  network  is  eo>*:'  to  the  value  of 
the  minimal  cut  set.  In  this  pa^er  the  value  of  a  cut  set  will  be 
referred  tc  as  its  capacity. 

C.  ENlMP.MiON  OF  CUT  SETS 

The  network  capacity  has  been  defined  to  be  equal  to  the  maximum 
flow  possible  In  the  network.  As  discussed,  this  maximum  flow  is  equal 
to  the  value  of  the  minimum  cut  set.  Therefore,  the  problem  of  mini¬ 
mizing  this  capacity  is  equivalent  to  minimizing  the  capac.ty  of  some 
cut  set.  Ii  is  obviou:  t  n«t  <">rtr«H  w.  ill  be  allocs  ted  tc  only  one  cut 
stt  since  If  this  were  not  the  case  all  aircraft  could  have  been  assigned 
to  the  cut  set  that  was  minimal  after  the  first  allocation  with  a 
resulting  decrease  in  network  capacity. 

The  complicating  factor  is  that  there  is  no  easy  way  to  find  out 
which  cut  set  should  be  selected  for  attack.  To  solve  the  problem  It 
is  necessary  to  have  some  means  of  identifying  cut  sets.  In  addition, 
it  is  desirable  to  be  able  to  identify  these  cut  sets  in  order  of 
Increasing  capacity  after  unlimited  interdiction  since  once  a  cut  set 
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;<  found  whose  capacity  after  unlimited  interdiction  is  greater  than  or 
equ*i  to  network  capacity  before  interdiction  no  more  cut  sets  need  be 
identified.  The  network  capacity  before  interdiction  represents  an 
upp^r  *ound  on  network  capacity.  Define  Sj  as  the  cut  set  with  the  i*h 
smallest  capacity  after  unlimited  interdiction.  The  set  of  S.  whose 
capacities  after  unlimited  interdiction  is  less  than  the  upper  bound  on 
network  capacity  will  be  denoted  by  S. 

The  method  by  which  cut  sets  are  identified  makes  use  of  the  topo¬ 
logical  dual  of  a  network.  Arcs  have  lengths  rather  than  capacities  in 
the  dual  network.  The  cut  sets  in  the  primal  network  have  a  one-to-one 
correspondence  with  the  loopless  paths  in  the  dual.  The  problem  of 
finding  the  shortest  path  from  the  dual  source  to  the  dual  sink  corre¬ 
sponds  to  the  primal  problem  of  finding  the  minimum  cut  set.  The  length 
of  the  dual  shortest  path  equals  the  primal  capacity. 

The  topological  dual  of  a  given  primal  network  is  constructed  as 
follows: 

(1)  Connect  the  source  and  the  sink  of  the  primal  with  an  artificial 
arc.  Call  the  result  the  modified  primal. 

(2)  Place  a  node  in  the  area  surrounding  the  modified  primal  (external 
face)  and  one  in  each  face  formed  by  the  arcs  of  the  modified  primal. 

Let  the  dual  source  be  the  node  in  the  external  face  and  the  dual  sink 
be  the  node  in  the  face  involving  the  artificial  arc. 

(3)  For  each  arc  In  the  primal  (except  the  artificial  arc)  construct 
a  dual  arc  that  Intersects  It  and  joins  the  two  nodes  in  the  faces 
adjacent  to  It. 

(4)  Assign  each  dual  arc  a  length  equal  to  the  capacity  of  the  primal 
arc  it  Intersects. 


Once  the  dual  network  has  been  developed,  the  shortest  path  through 
the  dual  before  interdiction  is  found.  This  path  is  determined  using 
the  upper  bounds  on  primal  capacities  as  lengths  of  arcs  in  the  dual. 

The  length  of  this  path  represents  network  capacity  before  interdiction. 
Any  shortest  path  algorithm  may  be  used  for  this  determination.  Dreyfus 
[4J  evaluated  several  of  these  algorithms  concluding  that  the  procedure 
developed  by  Dijkstra  is  the  most  efficient.  Next  the  lengths  of  the 
dual  arcs  are  changed  to  correspond  to  the  lower  bounds  on  primal  arc 
capacities.  The  lengths  of  the  dual  paths  now  represent  the  capacities 
of  the  corresponding  primal  cut  sets  after  unlimited  interdiction. 

Paths  with  loops  need  not  be  considered  since  they  correspond  to  primal 
cut  sets  that  either  include  more  arcs  than  necessary  to  sever  the 
network  or  contain  some  arc  more  than  once.  The  dual  paths  are  identi¬ 
fied  in  order  of  increasing  length  by  means  of  an  nth  shortest  path 
algorithm.  Clarke,  Krikorian  and  Rausen  [5]  developed  an  algorithm  for 
determining  the  n  best  loopless  paths,  but  it  is  difficult  to  apply. 
Pollack  [6]  in  an  unpublished  paper  presented  an  algorithm  which  succes¬ 
sively  develops  the  best  loopless  paths  using  extensions  of  shortest 
pat:;  algorithms.  This  procedure  is  less  complex  than  that  of  Clarke, 
Krikorian  and  Rausen  and  appears  to  be  more  efficient.  It  should  be 
noted  that  depending  on  the  number  of  elements  in  S  and  the  total 
number  of  paths  in  the  dual,  the  most  efficient  means  of  developing  S 
may  be  to  enumerate  all  paths  through  the  dual  and  then  compare  lengths. 
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III.  ANALYSIS  OF  THE  MODEL 


A.  MATHEMATICAL  FORMULATION 

The  problem,  as  outlined  previously,  is  to  find  that  allocation  of 
aircraft  to  an  airstrike  against  a  transportation  network  which  will 
minimize  the  capacity  of  that  network.  This  minimization  is  accomplished 
subject  to  a  constraint  on  aircraft  availability  and  the  consideration 
that  the  incremental  benefit  of  assigning  aircraft  must  exceed  the 
incremental  cost  of  that  assignment.  If  net  benefit  is  defined  to  be 
the  difference  between  the  total  benefit  derived  from  the  airstrike  and 
the  total  cost  of  aircraft  assignment  the  problem  may  be  restated  as 
follows:  maximize  net  benefit  subject  to  aircraft  availability. 

Let  K  represent  the  total  number  of  aircraft  available  for  assign¬ 
ment  to  the  airstrike  and  let  K*  be  the  number  of  aircraft  that  have 
been  assigned  to  the  airstrike.  Then  for  any  choice  of  K*  the  problem 
may  be  stated  mathematically  as 

min  [cut  set  capacity  after  optimal  interdiction] 

SiCS 

or 

min  [min  z  (1-.  +  w.  exp{-bi  .k.  })] 

S,cS  (1  »j)eS,  1J  ’•>  '•>  fJ 

subject  to  e  k* .  <  K* 

(i.ilcSj 

positive  integer  . 


( 

A 
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The  structure  of  this  problem  will  allow  the  development  of  an 
efficient  solution  procedure.  Note  that  with  respect  to  a  particular 
cut  set  the  objective  is  to  minimize  its  capacity.  Since  the  cut  set 
capacity  is  the  sum  of  functions  that  are  convex  in  k^,  this  capacity 
Is  a  convex  function  and  is  therefore  unimodal  with  respect  to  minimiza¬ 
tion.  The  overall  objective  function  is  the  minimum  of  a  set  of  convex 
functions  and  is  neither  concave  nor  convex.  This  together  with  the 
problem  of  not  knowing  which  cut  set  is  going  to  be  attacked  requires 
that  each  cut  set  in  S  be  the  subject  of  a  minimization  problem. 

For  a  particular  cut  set,  S-j,  the  problem  is 

min  1  (lii  +  w1i  expC-bi ik. .3) 

(i,j)eS.  1J  1J  10  1J 

subject  to  z  k.-.-  <  K 
Ci.jJeS,  ,J  - 

k^j  positive  integer  . 

The  term  z  1. .  is  constant  and  may  be  deleted  during  the  minimization 
<i,j)  1J 

and  then  added  back  to  give  the  solution  in  terms  of  capacity.  This 
problem  will  be  solved  by  means  of  dynamic  programming.  Each  arc  in 
the  cut  set  under  consideration  will  be  represented  by  a  stage  in  the 
dynamic  program. 

Let  the  number  of  arcs  be  n  and  resubscript  each  arc  (i,j)  and  its 
associated  parameters  in  any  order  with  the  single  subscript  i  running 
from  one  to  n.  The  decision  variable  for  stage  i  is  k^  and  the  return 
function  for  stage  i  is  given  by 

r1  *  wi  exp(-biki)  . 


14 


The  state  variable  for  stage  i  will  be  denoted  by  and  represents  the 
remaining  resource  availability  at  stage  i.  The  problem  may  be  restated 
as  finding  fn(xn)  where 


fn<xn> 


min  z  Mx. ) 
°ikni*n  1-’ 


subject  to  x^ =  x^  -  k..  . 

Fn(xn)  is  the  optimal  return  from  stages  n,n-l  ,...,1  given  xn  units  of 
resource.  The  above  problem  may  be  solved  by  dynamic  programming  since 
f n ( xn )  can  be  decomposed  into  a  series  of  single  variable  problems. 
Nemhauser  [7]  shows  that  problems  with  additive  stage  returns  may  always 
be  decomposed.  Therefore,  the  following  recursive  relationship  is  valid 


fn(xn)  =  min  Crn(kn)  + 

0<kn<_xn  ,  n>l 

and 

fi(Xi)  =  min  rWxi)  . 

°±kUxl 

The  value  of  xn_]  is  given  by 

xn-l  =  tn^xn»kn^  =  xn  "  kn 

where  tn  is  the  transformation  which  gives  the  relationship  between  the 
amount  of  resource  remaining  after  stage  n  given  that  xn  was  available 
before  stage  n  and  kn  was  utilized  at  stage  n. 

The  dynamic  program  is  solved  by  starting  at  stage  one  and  working 
to  stage  n  solving  a  series  of  single  variable  minimizations.  These 
minimizations  are  facilitated  by  the  convexity  of  the  individual  stage 
returns.  Nemhauser  [7]  provides  a  proof  of  the  fact  that  in  the 
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minimization  of  additive  stage  returns  the  convexity  of  each  stage 
return  ensures  that  f ^ )  is  a  convex  function  of  x^.  This  means  that 
each  single  variable  optimization  performed  in  the  dynamic  program  is 
of  a  unimodal  function  and  permits  the  use  of  Fibonacci  search  to  find 
the  optimal  values  of  the  decision  variables.  An  application  of  this 
technique  is  found  in  Ref.  7. 

After  the  optimal  allocation  of  aircraft  within  each  cut  set  in  S 
is  found  for  a  given  K*,  their  capacities  are  compared.  The  cut  set 
with  the  minimum  capacity  is  the  one  that  would  be  attacked  if  K*  air¬ 
craft  were  to  be  assigned  to  the  airs  trike.  The  capacity  of  this  minimal 
cut  set  is  by  definition  the  network  capacity  and  this  capacity  will  be  a 
strictly  decreasing  function  of  K*.  The  remaining  problem  is  to  deter¬ 
mine  how  many  aircraft  to  assign  to  the  strike  in  order  to  maximize 
net  benefit.  To  make  this  determination  it  is  necessary  to  know  the 
cost  of  allocating  aircraft  to  the  strike.  This  will  be  assumed  to  be 
a  constant  C  dollars  per  aircraft.  The  benefit  derived  from  network 
capacity  reduction  must  also  be  known.  It  will  be  assumed  to  be  a 
constant  D  dollars  per  unit  capacity  reduction. 

With  the  above  information  the  problem  of  determining  how  many 
aircraft  to  allocate  may  be  determined  by  comparing  the  incremental 
cost  of  assigning  aircraft  to  the  benefit  resulting  from  that  assign¬ 
ment.  To  make  this  comparison  it  is  necessary  to  define  the  benefit 
resulting  from  the  assignment  of  a  single  aircraft.  This  will  be  de¬ 
fined  as  the  product  of  the  benefit  per  unit  capacity  reduction  (D) 
and  the  amount  of  capacity  reduction  that  can  be  achieved  by  that 
aircraft. 
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The  amount  of  capacity  reduction  that  can  be  achieved  by  one 
additional  aircraft  is  a  function  of  the  number  already  assigned  and 
will  be  denoted  as  6 ( K* ) .  A  simple  decision  rule  is  to  assign  aircraft 
K*  *  1,2,...  until  a  point  is  reached  where  benefit  from  the  last 
aircraft  assigned  does  not  exceed  the  cost  of  assignment.  At  this  point 

D  *  6(K*)  <  C 
or 

S(K*)  £  C/D 

and  the  optimal  allocation  of  aircraft  is  K*-l  .  If  6 ( K*)  >  C/Q 

for  all  K*  the  optimal  allocation  is  K  under  this  rule.  There  would  be 

no  problems  with  this  decision  rule  if  6(K*)  were  a  non-increasing 

function  of  K*.  In  this  case  once  a  K*  was  found  such  that  6(K*)  C/D 

the  cost  of  any  further  assignment  of  aircraft  would  exceed  its  benefit. 

If  network  capacity  after  optimal  interdiction  were  determined  by 
only  one  cut  set  for  all  values  of  K*  then  s(K*)  would  be  non-increasing. 
This  is  not  the  case.  In  general  as  K*  ranges  from  0  to  K  different  cut 
sets  are  minimal  (see  Figure  1).  At  K*  =  0  the  cut  set  that  determines 
network  capacity  is  by  definition  the  one  that  is  minimal  before  any 
Interdiction  takes  place.  Unless  this  cut  set  is  also  minimal  after 
unlimited  interdiction,  at  some  point  another  cut  set  must  determine 
network  capacity.  This  crossover  may,  of  course,  occur  after  assignment 
of  all  available  aircraft.  These  points  where  a  change  in  the  constrain¬ 
ing  cut  set  occurs  represent  points  where  6(K*)  increases  with  respect  to 
K*.  Therefore,  there  is  no  guarantee  that  stopping  when  6(K*)  <_  C/D 
for  the  first  time  is  optimal.  If  at  some  point  after  further  assign¬ 
ment  of  aircraft  is  made  <5(K*)  again  exceeds  C/D,  it  may  be  that  further 
assignment  of  aircraft  would  have  resulted  in  benefits  outweighing  costs. 
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Capacity 


Figure  1.  Network  Capacity 


The  problem  of  determining  the  optimal  K*  will  be  handled  as  follows: 

(1)  Find  the  first  value  of  K*  for  which  5(K*)  <_  C/D.  Subtract  one 
aircraft  and  let  the  resulting  value  of  K*  be  K-j*.  If  K*  =  K  before 
K-|*  is  found  then  the  optimal  allocation  of  aircraft  is  K. 

(2)  Check  to  see  if  6 ( K* )  >  C/D  for  any  values  of  K*  >  Kj*.  If  not 
go  to  step  (4).  If  so  find  the  next  value  of  K*  for  which  6 (K*)  C/D. 
Let  this  number  minus  one  be  Kg*. 

(3)  Continue  in  this  manner  to  identify  the  K*  at  which  6(K*)  becomes 
<_  C/D  after  there  has  been  an  intervening  value  of  K*  such  that 

6 ( K* )  >  C/D.  Subtracting  one  aircraft  each  time,  label  the  resulting 

values  K3*,K4* . 1^*  If  s(K)  >  C/D  let  K,,*  =  K.  Let  Kq*  be 

defined  as  0. 

(4)  Starting  with  i  =  1  and  continuing  until  i  =  n,  check  whether  or 
not  the  cost  to  reach  K-j  from  K^.-j  is  exceeded  by  the  benefit.  If  it  is. 
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let  K*  =  Ki*  ,  increment  i  by  one,  and  go  to  the  beginning  of  step 
opt  1 

(4) .  If  it  is  not,  go  to  step  (5). 

(5)  Starting  with  1  =  1  and  continuing  until  1  =  n-i  check  whether  the 
cost  to  reach  from  K^_i  is  exceeded  by  the  benefit.  If  it  is,  let 
K*0pt  =  K*.j+j  ,  let  i  =  i+1+1  and  go  to  step  (4).  If  not,  increment  1 
by  one  and  go  to  the  beginning  of  step  (5). 

At  the  end  of  this  procedure  K*opt  will  be  the  optimal  number  of  air¬ 
craft  to  assign  to  the  airstrike  and  the  problem  will  be  solved. 

B.  STEPWISE  SOLUTION  PROCEDURE 

(1)  Formulate  the  topological  dual  of  the  transportation  network.  Find 
the  shortest  path  through  the  dual  before  interdiction.  This  represents 
an  upper  bound  on  network  capacity. 

(2)  Use  Pollack's  algorithm  [6]  to  identify  the  first,  second,  third, 
....  shortest  paths  through  the  dual  using  the  lower  bounds.  Continue 
identifying  paths  until  one  is  found  whose  length  exceeds  the  previously 
found  upper  bound  on  network  capacity.  Let  the  primal  cut  sets  corre¬ 
sponding  to  these  paths  be  denoted  as  set  S. 

(3)  For  each  cut  set  that  is  an  element  of  S,  use  dynamic  programming 
to  find  the  optimal  allocation  of  aircraft  and  the  resulting  capacity 
for  K*  equal  to  1,2,. ..  ,K. 

(4)  For  each  value  of  K*  find  the  network  capacity  by  taking  the  minimum 
of  the  capacities  of  the  elements  of  S. 

(5)  Construct  the  function  6{K*)  and  determine  K]*,  l<2*,...,Kn*. 

(6)  Using  the  procedure  previously  outlined  determine  which  K^*  is 
optimal . 


C.  SAMPLE  PROBLEM 

The  diagram  in  Figure  2  represents  a  hypothetical  transportation 
network.  The  three  numbers  associated  with  each  arc  are  bj j ,  1-jj, 
and  u^j. 


Figure  2.  An  Example  Network 


Figure  3  shows  how  the  topological  dual  of  the  transportation 
network  is  constructed. 


i 


Figure  3.  Construction  of  the  Dual 
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Figure  4  shows  the  topological  dual  after  the  data  for  each  arc 
has  been  transferred  from  the  primal.  In  the  dual  u^j  and  represent 
bounds  on  arc  length. 


Figure  4.  The  Topological  Dual 


To  simplify  notation,  paths  through  the  dual  will  be  designated  by  the 
nodes  over  which  they  pass.  The  shortest  path  through  the  dual  before 
Interdiction  is  1,2, 5,9  with  a  length  of  1395.  This  gives  an  upper  bound 
on  network  capacity.  Table  I  lists  all  loopless  paths  through  the  dual 
In  order  of  length  after  unlimited  Interdiction.  It  should  be  noted 
that  the  length  of  path  number  11,  the  11th  shortest  path  after  unlimited 
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TABLE  I.  DUAL  PATHS 


PATH 

NODES 

\ 

1.4. 7,6.9 

780 

2 

1.4, 7, 8.9 

890 

3 

14, 3, 2. 5.9 

960 

4 

1.2. .9 

1075 

5 

1.4.3.2.5.6.9 

1090 

6 

1.3.2, 5.9 

KrC 

7 

1.4.7.6.5.9 

1190 

8 

1.2.S.6.9 

1205 

9 

1.3.2.5.6.9 

1280 

10 

1,3.4, 7,6,9 

1290 

11 

1.3,4. 7. 0.9 

1400 

12 

1,4, 3,2,5, 6. 7. 8.9 

1600 

13 

1,2. 3,4, 7, 6.9 

1615 

14 

1.3,4, 7.6. 5. 9 

1700 

IS 

1.2. 5,6, 7. 8, 9 

1715 

16 

\2, 3.4, 7, 8. 9 

1725 

17 

1.3,2. 5,5, 7, 8.9 

1790 

18 

1.2, 3, 4, 7, 6. 5.9 

2025 
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Interdiction,  exceeds  the  upper  bound  on  network  capacity.  Therefore 
the  cut  sets  comprising  set  S  correspond  to  paths  1  through  10. 

It  is  assumed  for  purposes  of  this  example  that  there  are  100  air¬ 
craft  available  for  assignment  at  a  cost  of  30,000  dollars  for  each 
aircraft  assigned.  It  Is  further  assumed  that  the  benefit  derived  from 
a  reduction  of  one  ton  per  day  in  network  capacity  Is  7,500  dollars. 

The  dynamic  program  for  each  that  Is  constraining  Including  a 
sensitivity  analysis  on  K*  is  contained  In  the  computer  output.  A 
graph  of  the  resulting  network  capacity  is  given  by  Figure  5.  For  K* 

In  the  range  1  through  25  cut  set  4  determines  network  capacity,  for 
K*  In  the  range  26  through  57  cut  set  3  Is  constraining,  and  for  k*  from 
58  to  100  cut  set  1  Is  minimal. 

From  the  giver  values  of  C  and  D,  30,000  and  7,500  respectively, 
the  points  of  Interest  are  those  at  which  ;(*•)  becomes  C/0  »  4.  This 
occurs  for  the  first  time  when  k •  ■  42.  Therefore,  kj*  •  41.  At  k*  •  58 
4(k*)  again  exceeds  4  so  It  Is  necessary  to  search  for  another  point 
where  «(k#)  4.  This  next  occurs  at  k*  •  62  and  k^*  61.  Since 

5(K*)  does  not  exceed  4  for  any  k*  >  62,  k  •  ■  k^V 

It  Is  obvious  that  the  benefit  to  get  to  k^*  exceeded  the  cost 
since  kj*  was  the  first  point  at  which  the  allocation  of  another  aircraft 
did  not  produce  benefits  exceeding  costs.  However,  it  is  not  quite  as 
obvious  when  the  decision  Is  made  whether  or  not  to  allocate  k^* 
aircraft.  The  benefit  to  get  from  k^*  to  k^*  Is  equal  to  the  Incremental 
capacity  reduction  ■ultipiied  by  0. 

leneflt  •  62.23  *  7,500 
•  466,725. 
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This  benefit  Is  compared  to  the  cost  of  allocating  K2*  -  K^*  aircraft. 

Cost  =  20  X  30,000 
«=  600,000. 

Thus  the  benefit  is  outweighed  by  the  cost.  Since  there  is  no  allocation 
of  aircraft  greater  than  <2*  that  will  result  in  benefits  exceeding 
costs.  It  may  be  concluded  that  K*  =  41  represents  the  optimal  number 
of  aircraft  to  assign  to  the  airstrike.  At  this  level  of  interdiction 
the  network  capacity  will  be  1056.06  tons  per  day.  This  is  a  reduction 
of  338.94  tons  per  day  with  a  resulting  benefit  of  2,542,050  dollars. 

The  cost  of  this  reduction  is  1,230,000  dollars.  The  cut  set  that  will 
be  attacked  is  the  cut  set  corresponding  to  path  number  three  which 
contains  the  following  primal  arc*  ;  (4,7);  (4,6);  (2,6);  (1,6);  and 
(1,3).  These  arcs  correspond  to  dynamic  programming  stages  1,2, 3, 4, 
and  5  respectively.  Looking  at  the  dynamic  programming  stages  the 
optimal  allocation  of  aircraft  is:  k^j  =  9;  k4  6  =  6*  k2  6  =  7; 

*  10;  and  kj  3  *  9.  This  completes  the  solution. 
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IV.  DISCUSSION 


A.  PROPERTIES  OF  THE  SOLUTION  TECHNIQUE 

The  dynamic  programming  approach  taken  to  the  problem  guarantees 
that  the  solution  found  will  be  a  global  minimum  over  the  feasible 
region.  The  integer  constraints  pose  no  problem.  In  fact,  the  integer 
restriction  limits  the  number  of  values  the  decision  variables  may 
assume  and  allows  an  exact  solution  to  be  found.  Dynamic  programming 
also  provides  a  built  in  capability  for  sensitivity  analysis. 

The  convexity  of  the  damage  function  allowed  the  use  of  Fibonacci 
search  within  the  dynamic  program  resulting  in  a  tremendous  savings  in 
the  number  of  separate  calculations  made  in  each  dynamic  program.  With 
K  equal  to  100  the  reduction  was  on  the  order  of  10"^  times  the  number 
of  calculations  needed  for  exhaustive  search.  Larger  values  of  K  will 
produce  savings  of  an  even  greater  magnitude.  The  execution  time 
required  for  the  sample  problem  was  15.06  seconds  on  an  IBM  360/67. 
Utilizing  Fibonacci  search  it  was  found  that  the  increase  in  execution 
time  for  larger  values  of  K  was  approximately  linear.  Execution  time 
was  also  roughly  linear  with  respect  to  the  total  number  of  dynamic 
programming  stages  required  (45  in  the  sample  problem).  From  the  above 
observations  the  amount  of  computer  time  required  for  larger  problems 
may  be  predicted.  For  example,  a  problem  with  15  cut  sets  in  S  averaging 
6  arcs  per  cut  set  would  require  90  dynamic  programming  stages.  If  200 
aircraft  were  available  a  reasonable  estimate  would  be  that  this  problem 
would  take  approximately  4  times  as  long  to  solve  as  the  sample  problem. 

A  further  reduction  in  the  number  of  calculations  required  may  be 
achieved  with  a  coarse  grid.  Aircraft  can  be  allocated  in  packages  of 
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five  and  the  constraining  cut  sets  determined.  These  cut  sets  can  then 
have  aircraft  reallocated  one  at  a  time  and  the  optimal  solution  found 
as  before.  This  approach  can  not  guarantee  that  the  correct  constraining 
cut  sets  will  be  selected,  but  if  they  are  the  solution  will  be  optimal. 

Dynamic  programming  allows  some  generalizations  to  be  made  in  the 
problem.  To  begin  with,  since  additive  stage  returns  are  always 
decomposable,  the  technique  places  no  restrictions  on  the  damage  functions. 
The  negative  exponential  damage  function  used  in  this  paper  has  intuitive 
appeal  since  it  does  exhibit  diminishing  marginal  returns.  This  function 
also  contributes  to  ccmputational  efficiency  since  its  convexity  allowed 
the  use  of  Fibonacci  search.  However,  if  actual  interdiction  data 
suggests  damage  functions  of  another  form,  the  problem  can  still  be 
solved  with  somewhat  greater  expenditures  of  computer  time. 

Another  generalization  suggested  by  dynamic  programming  is  to  consider 
the  allocation  of  two  types  of  aircraft.  In  this  case  a  damage  function 
of  the  form 


*  'ij  tw1j  exp(-b,jkij  -  a,jh,j) 

might  be  assumed  with  k-jj,  1-jj,  and  w^j  defined  as  before,  h^  repre¬ 
senting  the  number  of  aircraft  of  the  second  type  assigned  to  arc  (i,j), 
and  denoting  the  vulnerability  parameter  corresponding  to  the 

second  type  of  aircraft.  Dynamic  programming  may  again  be  used  to  solve 
the  problem,  but  two  state  and  two  decision  variables  are  required. 

Although  the  new  damage  function  preserves  convexity,  in  this  case 
the  series  of  minimizations  is  of  functions  of  two  variables  and 
Fibonacci  search  is  not  applicable.  A  minimization  problem  was  run  for 
a  hypothetical  cut  set  containing  five  arcs.  The  execution  time  required 
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for  solution  was  5.63  seconds  when  10  aircraft  of  two  types  were 
available;  with  19  aircraft  of  each  type,  the  time  required  was  32.79 
seconds;  and  when  25  of  each  type  aircraft  were  available,  over  a 
minute  of  computer  time  was  used.  To  deal  with  even  relatively  small 
networks  the  computer  time  requirements  would  become  prohibitive  if  it 
were  necessary  to  consider  larger  aircraft  availabilities.  To  assign 
three  types  of  aircraft,  dynamic  programming  would  require  three  state 
and  three  decision  variables  and  the  technique  would  be  impractical 
even  for  small  problems. 

Another  application  of  the  dynamic  programming  approach  is  in  a 
modification  of  Nugent's  algorithm  [3].  This  modification  will  provide 
Integer  solutions.  Nugent  presented  a  method  of  finding  non-integer 
allocations  of  resources  that  would  minimize  network  capacity  subject 

to  i  k-n  <  K  and  k..  >  0  .  As  previously  discussed,  the 

(1,J)cS1  ij-  1J  " 

objective  function  in  this  problem  is  convex  with  respect  to  k^.  In 
Nugent's  formulation  the  feasible  region  defined  by  the  constraints  is 
also  convex.  Therefore,  for  any  particular  cut  set  the  problem  is  a 
convex  non-linear  program  and  Kuhn-Tucker  theory  provides  conditions 
that  are  both  necessary  and  sufficient  for  a  global  minimum.  Nugent 
solves  these  Kuhn-Tucker  conditions  and  using  an  upper  bounding  technique 
arrives  at  the  cut  set  that  will  be  minimal  after  optimum  interdiction. 

In  the  modification  the  set  S  and  the  upper  bound  on  network  capacity 
are  found  as  before.  The  Kuhn-Tucker  conditions  are  then  solved  to  find 
non-integer  constrained  solutions  that  minimize  cut  set  capacity  for 
each  element  of  S.  The  minimum  of  these  solutions  represents  the  optimal 
solution  without  integer  constraints.  When  integer  constraints  are  added 
this  minimum  represents  a  lower  bound  on  network  capacity.  The  cut  set 


with  the  minimal  non-integer  solution  is  deleted  from  S  and  becomes  the 
subject  of  a  dynamic  program  to  find  an  integer  solution.  If  this 
Integer  solution  is  less  than  the  non-integer  solutions  corresponding 
to  the  remaining  elements  of  S  it  is  optimal.  If  it  is  greater  than  some 
or  all  of  the  elements  of  S  it  represents  a  new,  smaller  upper  bound  on 
network,  capacity.  Any  elements  of  S  with  non-integer  solutions  greater 
than  this  new  upper  bound  are  deleted  from  S.  From  the  remaining  elements 
of  S  the  cut  set  with  the  smallest  non-integer  solution  is  selected  from 
S.  Again  dynamic  programming  used  to  find  a  new  integer  solution.  The 
new  integer  solution  is  compared  with  the  old  integer  solution  and  the 
minimum  is  called  the  current  integer  solution.  The  current  integer 
solution  is  then  compared  to  the  remaining  non-integer  solutions  and 
the  process  is  repeated.  This  iterative  procedure  is  continued  until 
either  S  is  the  null  set  or  until  the  current  integer  solution  is  less 
than  or  equal  the  non-integer  solutions  corresponding  to  all  of  the 
remaining  elements  of  S.  In  either  case  the  current  integer  solution 
represents  the  optimal  solution  to  the  integer  constrained  problem. 

In  general,  if  the  number  of  aircraft  to  be  allocated  to  the  air- 
strike  is  known,  this  modification  is  more  efficient  that  using  dynamic 
programming  on  every  element  of  S  to  solve  the  minimization  problem. 

In  solving  the  example  problem  from  Nugent's  paper  it  was  necessary  to 
run  only  one  dynamic  program  and  with  exponential  damage  functions  the 
Kuhn-Tucker  conditions  are  easy  to  solve  relative  to  solving  a  dynamic 
program.  However,  this  modification  does  not  lend  itself  to  the 
sensitivity  analysis  on  K  that  is  necessary  when  the  number  of  aircraft 
to  be  assigned  to  the  strike  is  taken  to  be  a  decision  variable. 
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As  already  mentioned,  there  are  limitations  on  the  technique 
presented.  One  difficulty  that  has  not  yet  been  discussed  is  in  the 
measurement  of  the  costs  and  benefits  of  aircraft  assignment.  In  this 
paper  the  problem  was  ignored  and  constant  dollar  values  of  C  and  D 
were  selected  arbitrarily.  This  problem  is  important  since  the  selection 
of  C  and  D  determines  how  many  aircraft  will  be  assigned  to  the  strike. 

i 

If  D  had  been  taken  to  be  D00  dollars  per  ton  of  flow  reduced  vice 
7,500  and  the  rest  of  the  pr  'em  remained  unchanged,  the  decision  would 
have  been  made  to  allocate  77  aircraft  in  a  strike  against  cut  set  one 
resulting  in  a  network  capacity  of  938.19  tons  per  day.  On  the  other 
hand,  if  D  was  less  than  1,519  dollars  per  ton  of  flow  reduced  the 
solution  would  be  to  make  no  attack  against  the  network. 

B.  RECOMMENDATIONS  FOR  FURTHER  STUDY 

The  possibility  of  deriving  damage  functions  from  actual  interdiction 
data  was  mentioned  earlier.  If  the  method  of  this  paper  were  to  be  put 
to  use  in  solving  a  real-world  interdiction  problem  some  verification 
of  the  damage  function  would  be  essential.  However,  due  to  the  sensi¬ 
tivity  of  the  solution  to  both  costs  and  benefits,  the  measurement 
problem  associated  with  costs  and  benefits  should  receive  at  least  as 
much  attention  as  the  damage  function. 

Another  possibility  for  further  study  would  be  the  utilization  of 
the  model  described  in  this  paper  to  represent  real-world  problems  other 
than  aircraft  interdiction.  One  obvious  example  might  be  the  problem  of 
allocating  resources  to  the  improvement  of  a  highway  system.  In  this 
example  It  would  probably  be  relatively  easy  to  get  data  from  which 
to  derive  improvement  functions,  but  the  measurement  of  costs  and  benefits 
would  be  as  difficult  as  before. 
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The  model  presented  could  be  refined  by  assigning  different  values 
to  capacity  reduction  in  the  various  arcs  of  the  network.  The  objective 
then  would  be  to  minimize  the  maximum  value  of  flow  possible  in  the 
network  rather  than  to  minimize  network  capacity.  The  solution  technique 
presented  could  still  be  used.  A  further  refinement  might  be  to  consider 
not  only  arc  vulnerability,  but  also  the  repair  capability  of  the 
opponent.  This  would  require  capacity  reduction  to  be  taken  as  a 
function  of  time  as  well  as  aircraft  allocation  and  would  make  the 
analysis  of  the  model  more  difficult.  Many  other  refinements  could  be 
made  in  order  to  make  the  model  more  representative  of  the  real  world, 
but  in  general  the  increased  realism  gained  would  be  at  the  expense  of 
Increased  computational  effort. 


V.  SUMMARY 


A  solution  procedure  has  been  developed  for  the  problem  of 
determining  the  optimal  allocation  of  aircraft  in  planning  an  airs trike 
against  a  transportation  network.  The  damage  function  for  arcs  in  the 
network  is  assumed  to  have  a  negative  exponential  form.  To  make  use 
of  the  procedure  it  is  necessary  to  have  available  the  following  infor¬ 
mation:  the  upper  and  lower  bounds  on  the  capacity  of  each  arc,  the 
vulnerability  parameter  for  each  arc,  the  number  of  aircraft  available 
for  assignment  to  the  airstrike,  the  cost  of  assigning  an  aircraft  to 
the  strike,  and  the  benefit  resulting  from  network  capacity  reduction. 

In  the  solution  procedure  every  cut  set  that  is  designated  a 
candidate  for  attack  is  the  subject  of  a  dynamic  program.  A  sensitivity 
analysis  is  performed  on  the  number  of  aircraft  to  be  assigned  and  this 
gives  the  network  capacity  after  optimal  interdiction  as  a  function  of 
the  number  of  aircraft  assigned  to  the  strike.  A  cost  benefit  analysis 
Is  then  made  to  determine  the  largest  number  of  aircraft  that  can  be 
assigned  before  costs  of  further  allocation  begin  to  outweigh  the 
benefits  resulting  from  that  allocation. 

At  the  end  of  the  procedure  the  solution  consists  of  the  following: 
the  number  of  aircraft  to  assign  to  the  airstrike,  the  cut  set  that  will 
be  attacked,  the  number  of  aircraft  to  allocate  to  each  arc  of  the  cut 
set  chosen,  and  the  capacity  of  the  network  after  this  assignment  of 
aircraft. 
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COMPUTER  OUTPUT 


DP  SOLNS  CUT  SET  1 

ACFT  AVAIL  CAPACITY  CHANGE  IN  CAPACITY 


0 

1500.0000 

0.0 

1 

1472.2615 

-27.7385 

2 

1450.4417 

-21.8198 

3 

1430.5337 

-19.9080 

4 

1413.3696 

-17.1641 

5 

1 397. 6652 

-15.5044 

6 
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-13.5017 

7 

1372.2888 

-12.0747 

8 
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-11.7632 

9 
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-11.3020 

10 
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-10.8589 

11 
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-10.6206 

12 
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13 
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14 
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15 
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16 
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17 
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18 
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19 
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-8.8906 

20 

1242.2502 

-8.8259 

21 
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22 

1225.3132 

-8.3953 

23 

1216.9587 

-8.3545 

24 

1208.7517 

-8.2070 

25 

1200.7656 

-7.9861 

26 

1192.8806 

-7.8850 

27 

1165.2842 

-7.5964 

28 

1177.7083 

-7.5759 

29 

1170.3845 

-7.3237 

30 

1163.1057 

-7.2788 

31 

1155.8796 

-7.2261 

32 

1148.8862 

-6.9934 

33 

1142.0127 

-6.8735 

34 

1135.2935 

-6.7192 

35 

1128.7214 

-6.5720 

36 

1122.1831 

-6.5383 

37 

1115.7273 

-6.4558 

38 

1109.5078 

-6.2195 

39 

1103.3052 

-6.2026 

40 

1 C  57.  3457 

-5.9595 

41 

1091.4297 

-5.9160 

42 

1085.7039 

-5.7258 

A3 

1C  80. 0000 

-5.7039 

44 
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-5.6274 

45 

1068.8713 

-5.5012 

46 
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47 
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48 
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49 
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92 

93 

94 
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1033.1694 
1028.4817 
1023.8743 
1019.3701 
1014.9280 
1010.5452 
1006.2178 
1002.0488 
997.8911 
993.8245 
989.  8296 
985.8640 
982.0259 
978.2537 
974.5659 
970.9778 
967.4346 
963.9751 
960.5618 
957.1577 
953.8872 
950.6404 
947.4414 
944.2991 
941.2104 
938.1914 
935.2534 
932.3528 
929.5581 
926.7710 
924.0769 
921.39  9  2 
918.7407 
916.1680 
913.6394 
911.1230 
908.6511 
906.2458 
903.8708 
901.5830 
899. 3010 
897.1086 
894.9324 
692.  8259 
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888.6572 
686.6335 
884.6541 
882.6848 


CHANGE  IN  CAPACITY 

-5.0784 
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-4.3828 

-4.3274 
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-1.9794 

-1.9693 
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STAGE  NUMBER  A 
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10 

A 

11 

13 

Q 

14 

16 

l 

17 
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2 

20 

22 

A 

23 

25 

5 

26 

28 

6 

29 

31 

7 

32 

3A 

8 

35 

37 

9 

38 

AO 

LG 

41 

A3 

1 

AA 

A6 

13 

A7 

A9 

A 

50 

52 

15 

53 

55 

6 

56 

58 

7 

59 

61 

18 

62 

6A 

19 

65 

67 

21 

68 

70 

22 

71 

73 

23 

74 

76 

2A 

77 

79 

25 

80 

82 

26 

83 

85 

27 
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30 

92 

9A 
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95 

97 
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98 
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2A 
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DP  SOLNS  CUT  SET  3 


ACFT  AVAIL  CAPACITY  CHANGE  IN  CAPACITY 


0 

1535.0000 

0.0 

1 

1497.0632 

-37.9368 

2 

1469. 3250 

-27.7383 

3 

1443.3813 

-25.9436 

4 

1421.5618 

-21.8196 

5 

1401.7942 

-19.7676 

6 

1384.0525 

-17.7417 

7 

1366.6409 

-17.4116 

8 

1349.4768 

-17.1641 

9 

1332.9653 

-16.5115 

10 

1317.  1670 

-15.7983 

11 

1302.1807 

-14.9863 

12 

1288.3894 

-13.7913 

13 

1274. 8879 

-13.5015 

14 

1261.9890 

-12.8989 

15 

1249.3105 

-12.6785 

16 

1237.1775 

-12.1331 

17 

1225.6582 

-11.5193 

18 

1214.5559 

-11.1021 

19 

i; 

203.9353 

-10.6208 

20 

11 

L93.7605 

-10.1748 

21 

I  1 

84.  1384 

-9.6219 

22 

V. 

L74.5828 

-9.5556 

23 

11 

L66.2283 

-8.3547 

24 

11 

L5  7. 9309 

-8.2972 

25 

11 

L49.7063 

-8. 2246 

26 

11 

41 .5410 

-8.1655 

27 

11 

L33.5042 

-8.0369 

28 

i: 

L26.4250 

-7.0790 

29 

1119.7122 

-6.7130 

30 

1113.1401 

-6. 5720 

31 

1106.5872 

-6.  5529 

32 

1100.4941 

-6.0930 

33 

1094.8201 

-5.6741 

34 

1089.2129 

-5.6072 

35 

1083.9541 

-5.2589 

36 

1078.7097 

-5.2442 

37 

1073.5400 

-5.1697 

38 

1068. 8567 

-4.6835 

39 

1064.34:8 

-4.5138 

40 

1C 

560.1226 

-4.2203 

41 

1< 

556.0559 

-4.0666 

42 

1052.1440 

-3.9120 

43 

1048.2588 

-3.8850 

44 

1044.3787 

-3.8803 

45 

1040.°917 

-3.3869 

46 

1037.6477 

-3.3439 

47 

1034.3801 

-3.2675 

48 

1031.1814 

-3.1989 

49 

1028.3032 

-2.8781 

40 


ACFT  AVAIL 


CAPACITY 


CHANGE  IN  CAPACITY 


50 

1025.5740 

-2.7293 

51 

1022.8560 

-2.7181 

52 

1020.2021 

-2.6536 

53 

1017.6858 

-2.5164 

54 

1015.2087 

-2.4772 

55 

1012.9290 

-2.2797 

56 

1010.7476 

-2.1813 

57 

1CC8.6155 

-2.1321 

58 

1006.6360 

-1.9794 

59 

1004.7319 

-1.9042 

60 

1002.8967 

-1.8352 

61 

1001.0820 

-1.8147 

62 

999.3315 

-1.7505 

63 

997.7410 

-1.5905 

64 

996.1616 

-1.5795 

65 

994.6045 

-1.5571 

66 

993.  1997 

-1.4048 

67 

991.8401 

-1.3595 

68 

990.5115 

-1.3285 

69 

989.2705 

-1.2410 

70 

988.0457 

-1.2249 

71 

586.8755 

-1.1702 

72 

985.7483 

-1.1274 

73 

984.6384 

-1.1096 

74 

983.6313 

-1.0072 

75 

982.6680 

-0.9635 

76 

581.7410 

-0.9269 

77 

980.8362 

-0.9048 

78 

979.9695 

-0.8669 

79 

579.1208 

-0.8487 

80 

978.3464 

-0.7742 

81 

577.5886 

-0.7579 

82 

576.8425 

-0.7461 

83 

976.1165 

-0.7261 

84 

975.4697 

-0.6466 

85 

974.8276 

-0.6422 

86 

974.2314 

-0.5962 

87 

573.6487 

-0.5827 

88 

573.0684 

-0.5804 

89 

972.5156 

-0.5527 

90 

571.  9753 

-0.5401 

91 

571.4998 

-0.4757 

92 

571.0308 

-0.4690 

93 

570.5630 

—0*4676 

94 

970.1121 

-0.4511 

95 

969.7024 

-0.4095 

96 

969.3057 

-0.3969 

97 

968.9287 

-0.3768 

98 

568.5535 

-0.3753 

99 

568.1846 

-0.3689 

100 

967.8320 

-0.3524 
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95 

98 


0 

0 

1 

2 

3 

5 

6 
7 
9 
0 
1 

12 

13 

4 
6 

17 

18 
19 
21 
22 

23 

24 

25 

26 
28 

29 

30 

31 

33 

34 

35 

36 

37 
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STAGE  NUMBER  A 


1  STATE  I  DEC  I  STATE  IDEC  1ST  ATE 


IOEC 
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STAGE  NUMBER  5 


|  STATE  IDEC  1ST  AT  6 


IOEC  1  STATE  1DEC 


0 

3 

6 

9 

12 

15 

18 

21 

24 

27 

30 

33 

36 

39 

42 

45 

48 

51 

54 

57 

60 

63 

66 

69 

72 

75 

78 

81 

84 

87 

90 

93 

96 

99 


0 

0 

1 

2 

3 

3 

4 

5 

5 

6 
7 

7 

8 
9 

10 

10 

11 

12 

12 

13 

14 

15 

15 

16 
16 

17 

18 

19 

20 
20 
21 
21 
22 
23 


1 

4 

7 

10 

13 

16 

19 

22 

25 

28 

31 

34 

37 

40 

43 

46 

49 

52 

55 

58 

61 

64 

67 

70 

73 

76 

79 

82 

85 

88 

91 

94 

97 

100 


0 

0 

1 

2 

3 

3 

4 

5 

5 

6 

7 

8 
8 
9 

10 

10 

11 

12 

13 

3 

4 

5 

5 

6 

17 

18 
18 
IS 
20 
20 
21 
22 
23 
23 


2 

5 

8 

11 

14 

17 

20 

2? 

26 

29 

32 

35 

38 

41 

44 

47 

50 

53 

56 

59 

62 

65 

68 

71 

74 

77 

80 

83 

86 

89 

92 

95 

98 


0 

1 

1 

2 

3 

4 

4 

5 
5 
7 

7 

8 
9 
9 

10 

11 

12 

12 

13 

14 

14 

15 

16 
16 

17 

18 

19 

\l 

20 
21 
22 
23 
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OP  SOLNS  CUT  SET  4 


ACFT  AVAIL  CAPACITY  CHANGE  IN  CAPACITY 


0 

1395.0000 

0.0 

1 

1375.2324 

-19.7676 

2 

1357. 8208 

-17.4116 

3 

1341.3096 

-16.5112 

4 

1326.3232 

-14.9861 

5 

1312.5320 

-13.7914 

6 

1298.9368 

-13.5952 

7 

1286.0381 

-  2.8988 

8 

1274.5186 

-11.5195 

9 

1263.3877 

-11.1208 

10 

1252.2856 

-11.1021 

11 

1242.6638 

-9.6219 

12 

1233.1082 

-9.5556 

13 

1223.9949 

-9.1131 

14 

1215.7703 

-8.2246 

15 

1207.7334 

-8.0369 

16 

1200.2722 

-7.4612 

17 

1153.1934 

-7.0790 

18 

1186.4802 

-6.7130 

19 

1 EO. 3716 

-6.1087 

20 

1174.2786 

-6.0930 

21 

1168.6714 

-5.6072 

22 

1163.4272 

5.2442 

23 

1158.4258 

-5.0014 

24 

1153.7424 

-4.6835 

25 

1149.2285 

-4.5138 

26 

1145.1338 

-4.0948 

27 

1141.2219 

-3.9120 

28 

1137.3367 

-3.8850 

29 

1133.9841 

-3.3525 

30 

1130.6404 

-3.3439 

31 

1127.3728 

-3.2675 

32 

1124.4946 

-2.8781 

33 

1121.7500 

-2.7448 

34 

1119.0205 

-2.7293 

35 

1116.5435 

-2.4772 

36 

1114.2637 

-2.2797 

37 

1112.0164 

-2.2473 

38 

1109.8843 

-2.1322 

39 

11C7.9800 

-1.9042 

40 

1106.1401 

-1.8399 

41 

1104.3049 

-1.8352 

42 

1102.7146 

-1.5905 

43 

1101.1350 

-1.5795 

44 

1099.6287 

-1.5064 

45 

1098.2690 

-1.3595 

46 

1096.9407 

-1.3285 

47 

1055.7073 

-1.2333 

48 

1094.5371 

-1.1702 

49 

1093.4275 

-1.1096 

47 


r 


ACFT  AVAIL  CAPACITY  CHANGE  IN  CAPACITY 


50 

1052*4177 

-1.0098 

51 

1091.4106 

-1.0072 

52 

1090.4836 

-0.9269 

53 

1089.6169 

-0.8669 

54 

1088.7900 

-0.8267 

55 

1C  88. 0159 

-0.7742 

56 

1087.2698 

-0.7461 

57 

1086.5930 

-0.6769 

58 

1065.9463 

-0.6466 

59 

1085.3042 

-0.6422 

60 

1084.7500 

-0.5542 

61 

1C84.1973 

-0.5527 

62 

1083.6570 

-0.5401 

63 

1083.1814 

-0.4757 

64 

1082.7275 

-0.4537 

65 

1082.2764 

-0.4511 

66 

1081.8669 

-0.409  5 

67 

1081.4902 

-0.3768 

68 

1081.1187 

-0.3715 

69 

1 C  80. 766 1 

-0.3524 

70 

1080.4514 

-0.3148 

71 

1080.1472 

-0.3041 

72 

1C75.8440 

-0.3033 

73 

1079.5811 

-0.2629 

74 

1079.3201 

-0.2611 

75 

1079.0710 

-0.2490 

76 

1078.8462 

-0.2247 

77 

1078.6267 

-0.2196 

78 

1078.4229 

-0.2039 

79 

1078.2292 

-0.1934 

80 

1C78.0455 

-0.1834 

81 

1077.8789 

-0.1669 

82 

1077.7126 

-0.1665 

83 

1C77.5593 

-0.1532 

84 

1077.4160 

-0.1433 

85 

1077.2793 

-0.1367 

86 

1C77. 1514 

-0.1280 

87 

1077. 0281 

-0.1233 

88 

1076.9163 

-0.1119 

89 

1076.8093 

-0.1069 

90 

1076.7031 

-0.1062 

91 

1C76.6116 

-0.0916 

92 

1076.5203 

-0.0914 

93 

1076.4309 

-0.0893 

94 

1076.3523 

-0.0786 

95 

1076.2773 

-0.0750 

96 

1076.2026 

-0.0746 

97 

1076. 1350 

-0.0677 

98 

1076.0728 

-0.0623 

99 

1C76.0112 

-0.0614 

100 

1075.9531 

-0.0583 
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STAGE  NUMBER 


1 


I  STATE 


I  DEC 

1ST  ATE 

IDEC 

ISTATE 

IDEC 

0 

1 

1 

2 

2 

3 

4 

4 

5 

5 

6 

7 

7 

8 

8 

9 

10 

10 

11 

11 

12 

13 

13 

14 

14 

15 

16 

16 

17 

17 

18 

19 

19 

20 

20 

21 

22 

22 

23 

23 

24 

25 

25 

26 

26 

27 

28 

28 

29 

29 

30 

31 

31 

32 

32 

33 

34 

34 

35 

35 

36 

37 

37 

38 

3e 

39 

40 

40 

41 

41 

42 

43 

43 

44 

44 

45 

46 

46 

47 

47 

48 

49 

49 

50 

50 

51 

52 

52 

53 

53 

54 

55 

55 

56 

56 

57 

58 

58 

59 

59 

60 

61 

61 

62 

62 

63 

64 

64 

65 

65 

66 

67 

67 

68 

68 

69 

70 

70 

71 

71 

72 

73 

73 

74 

74 

75 

76 

76 

77 

77 

78 

79 

79 

80 

80 

81 

82 

82 

83 

83 

84 

85 

85 

86 

86 

67 

88 

88 

89 

89 

90 

91 

91 

92 

92 

93 

94 

94 

95 

95 

56 

97 

97 

98 

98 

99 

100 

100 
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STAGE  NUMBER  2 


ISTATE 


I  DEC  ISTATE 


I  DEC  ISTATE 


I  DEC 


50 


STAGE  NUMBER  3 


1ST  ATE  ICEC  ISTATE  IDEC  I  STATE 


I  DEC 


COMPUTER  PROGRAM 


THIS  PROGRAM  IS  DESIGNED  TO  FIND  THE  MINI MU*4  OF  A 
SEQUENCE  OF  FUNCTIONS.  EACH  FUNCTION  IN  THE  SEQUENCE 
IS  A  SUM  OP  NEGATIVE  EXPONENTIAL  FUNCTIONS  OF  Tw  = 

F CRM  BLO*W»ExP( !D£C) •  IOEC  REPRESENTS  THE  AMOUNT 
OF  RESOURCE  ALLOCATED  TO  REDUCE  THE  VALUE  CF  A 
PARTICULAR  NEGaTIV*  EXPONENTIAL  AND  IS  THE  DECISION 
VARIABLE.  THE  MINIMIZATION  tS  SUBJECT  TO  CONSTRAINT 
THAT  THe  SUM  op  THE  IDEC'S  FOR  EACH  FUNCTION  IN  THE 
SEQUENCE  IS  LESS  Than  OP  EOUAL  TO  k  whEPE  k 
REPRESENTS  Th|  TOTAL  RESOURCE  AVAUABlLTY.  DYNAMIC 
PROGRAMMING  IS  THE  SOLUTION  TECHNIQUE  USED. 

THE  INPUT  PARAMETERS  FOR  THIS  PROGRAM  ARE  AS  FOLLOWS  - 

NCUT -NUMBER  OF  FUNCTIONS  TO  P*  MINJMIZED. 

K  -  MAXIMUM  amount  OF  RESOURCE  AVAILABLE. 

UPBNO  -  PR£  Dc  TE°  M  IN^C  U°P  F R  BOUND. 

n  -  number  o*  Exponentials  in  particular  fun.tion. 

B.W.PLO  -  CONSTANTS  ASSOC  WITH  EACH  EXPONENTIAL. 


THE  OUTPUT  FROM  THIS  PROGRAM  IS  AS  FOLLOWS  - 

*A  EACH  FUNCTION  N  DYNAMIC  PROGRAMMING  STAGES 
REQUIRED.  C0R  PACH  STAGE  TK" 

MAL  VALUE  OF  TH t  DECISION  VARIABLE  IS  PRINTED 

For  fach  possible  value  of  the  state  variable.! 


UU  I  r 

i£iR 

OPTI 


I  STATE  -  VALUE  op  STATE 
IOEC  -  OPTIMAL  VALUE  FO 


VAR! Adi c. 

R  DEC!  SI  ON  VARIABLE. 


1 AT  THE  END 

printed  for 


JF  T*E  STAGES  THE  SOLUTION  IS 
■ACM  FUNCTION! 


AIRCRAFT  AVAIL  -  AMOUNT  OF  P'SPURCE  AVAILABLE. 
CAPAjjTY  r  SOIJTION  TO  MJ  *.’I  VI  ZA  TJ  ON  PRC6Lfv- 


DIM 

1986 


IN  CF.-  ENTAL  I MPROV  L  “tNT  IN  SOLUlfoN, 

MFNSION  FCN(?COI.FC  "E  w!  200)  .CAP  (  2001 
MANSION  !PL  A*"-  I20P  t  <  PL  AN  t  (  ZOO  I ,  DcLCAP  I  200 ) 
MeNSION  IflBNJCOi 

TA  JF  IBNO/  1.  2.4.  7  •  12,20.  33.54,88.143,232,376,609, 
6, 1596, 25?3, 4190. 6764, 10945, 1/710/ 


•  PLANES  AVAIL 


READ  4  PATHS  IN  NETWORK, 

4  LUB  ON  C4PAC1  TV 

PE  AO  (5 ,40  I  NCUT  ,  v  .UPBNO 
40  FORMA  T  (  2  if  3)  ,F9.*| 

DO  1200  IP  ATM* 1. NCUT 

DATA  FCN/2003.0/,  FCNNEw/  200*0.0/.  CAP/20  0*0.0/ 
DATA  I  PL  AN?  /  ZOO*  0  /  ,  «PL  ANE  /200*  0  / 

DATA  DELCAP/Z0U*0.0/ 

READ  4  ARCS  IN  CUT  SET 

PEA0I5.5O)  N 
50  FORMAT!  I  2) 

M-KM 

READ  IN  ARC  PARAMETERS  (1ST  ARC) 


PE  AO  1 5  .601  B.W.BLO 
60  FORMAT! 3(F 10.6) ) 
8LOSUM.BLO 

OP  STAGE  1 


BciSWii  i 
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FCN(  I)  =W*EXP(ARG) 

I  PLANE  ( 1  >=I-1 
KPLANE ( I )  = IPLANE( 1 1 
100  CONTINUE 
ISTAGE= 1 

WRITE (6, 1000 )  I  STAGE* (KPLANE( I ) 1 1  PLANE  (  I )  * 
OP  STAGES  2  THR U  N 
00  900  I STAGE*2  *  N 

READ  IN  ARC  PARAMETERS  (REMAINING  ARCS) 

READ (5*60)  B.W.BLO 
BLOSUM=B  LOSUM*  8LC 
IF(BLOSUM.GE.UPBNO)  GO  TO  2550 

SET  STATE  VARIABLE  £  RUN  FIB  SEARCH 

00  800  J=2.M 
IX«J-1 

DO  200  NO* 1,20 

IFUFIBNO(NO).GE.J)  GO  TO  300 
200  CONTINUE 
300  NCM*N0-1 

1b*IFIBN0(N0)-1 
DO  500  ITER  =  1, NQM 
N01*N0-ITER 
N02=NQ- 1  TER- 1 

IF(IA4IFIBNC(N01).G5.J )  GO  TO  400 
ARG1*-B*  (I  4+! p I  BNO ( N01  )  ) 
IT1=IX-IA-IFIBN0(N01)41 
!f(NC2.=Q.O)  GO  TC  350 
ARG2*-BMIA4lFIBN0<NC2  )) 
IT2=IX-IA-IFIBNO(N02)  +  1 
GO  TO  360 
350  ARG2*-B4 I A 
IT2=IX-!A*1 

360  Fl*W*cXP(APGl)4FCN< IT1 ) 

F2*W>  f  XP(ARG2)4FCNUT2  ) 

IF(P2.lT.F1)  GO  TO  400 
IA*IFIBNC(N02)4l4l  A 
GO  TO  500 

400  IB*IFIBN0(N01 J-l+IA 
500  CONTINUE 

rCNNEW  (J)*AMlNl(F2fFl) 

I PLANS ( J ) *  I  A 

I F (  I  STAGE. EQ.N)  CA P (  J)  *FCNNEW ( J) +  BLOSUM 
IF(J$TAGE.  EQ.N.AN0.J.GT.2) 

1DELCAP( J)=FCNMcW( J)-PCNNEW(IX) 

800  CONTINUE 

FCN(1 )*FCN(1 )♦« 

DO  850  1  *2  , M 
FCN(  I  )  =  FCNNEW(  I ) 

850  CONTINUE 

WRITE (6, 1000)  I  STAGE, (KPLANE (I ) , I  PLANE  (I  >  , 
900  CONTINUE 

CAP(1  )  =FCN( 1 l  +  BLOSUM 
OELCAP (  2  )  *FCN(  2  )  -PC  N (1 ) 


ENO  OP  STAGES 


1000  FORMAT ('l*//////// 
1*  •  ,14X,3( • 1ST  ATE* 
14X ,  3(  I6,4X,I4 
L4X. 3( I6.4X, 14 
:4X  ,3  (  16, 4X,  14 
14  X  ,  3(  I 6,4X* I  4 
14X  ,  3 (  16,  4X,  14 
L4  X  ,  3  (  16, 4X,  14 


,2X1  /• 
, 2X ) /  • 
.2X1/* 
,2X)  /• 
, 2X) /* 
,2 X )  / • 
•  14X,  3(  I6.4X.I4.2X)  /• 


•  • »29X»  • 
,4X.  MOEC 

- I  I  l  ■ 

• 
t 
t 
f 
f 

9 


STAGE  NUMBER  * . 
• ,2X)/// 

14X , 3 ( 16 ,4X ,  I  4, 
14X, 3( I  6, 4X, I  4 , 
14X,  3 (  16,  4X,  14, 
14X,3(I6,4X.I4, 
1 4X , 3(  I  6, 4X, 1  4 , 
14X,3( 16, 4X, 14, 
14X.3(  I  6.4X.I4  . 


1=1, M) 


1=1, M) 


12///// 

2X  )/ 

2  X )  / 

2  X  )  / 
2X1/ 

2  X )  / 

2X  )/ 

2X)  / 
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,14X, 3( 16, 4X, 14, 2X1/ 
,14X,3(I6,4X,I4,2X)/ 
.14X.  3U6.4X.I4.2X)  / 
t 14X»  3( 16, 4X, I4»  2X )/ 
,14X,3(I6,4X»I4»2X)/ 
,  14X, 3‘ I 6,4X,I4,2X)/ 
, 14X , 3  ( 16, 4X »  14, 2X  )/ 
» 1 4  X  •  3  (  I6.4X.I4.2X)  / 
,14X,3< 16 . 4X , I  4 , 2X )  / 
,14X,3(I6,4X,I4,2X)/ 
.14X.3(I 6.4X.I4.2X)  ) 


1050 


»27X ,  •  DP 


•,5X, 'CAPAC I TY* , 


SCLNS  CUT 
5X, 


1075 

1100 

1125 


1150 

1200 

2550 

2575 

2600 


,14X, 3(16, 4X, 14,  2X )  /  ‘ 
,14X,3(I6,4X,I4,2X)/' 
.14X.3(I6.4X,I4,2X)/, 
,14X,3(I6,4X,I4,2X)/» 
,14X,3(I6,4X,I4,2X>/' 

*  14X ,3(  16, 4X,  14,  2X ) /  ' 
,14X,3(I6,4X,I4,2X)/' 

,  14X.  3(  I6.4X.I4. 2X)  /• 

,  14X  ,3(I6,4X,I4,2X)/' 
,14X,3(I6,4X,I4,2X)/' 

. 14X, 3( 16, 4X.I4.2X)/* 

WRITE(6 ,1050  )  IPATH 
FORMAT ( ‘l* ////////•  • 

L  •  '. 12X.  'ACFT  AVAIL 
['CHANGE  IN  CAPACITY'//) 

M12=M/2 
M121=M12+1 
DO  1100  1=1,  Ml  2 

WRITE ( 6. 1075)  KPLANE( I ) ,CAP( I ) .DELCAP(I) 
FORMAT (•  • ,10X,I10,3X,F10.4, 13X,F10.4) 
CONTINUE 
WRITE (6,1125) 

FORMAT  Cl'//////// 

['  ' ,  12  X ,  'ACF  T  AVAIL' ,5X, 'CAPACITY*  ,5X, 
['CHANGE  IN  CAPACITY'//) 

DO  115C  I =M121 , M 

WRITE( 6. 1075)  KPLANE(I) ,CAP( I ) .DELCAP(I) 
CONTINUE 
CONTINUE 
GO  TO  2600 
WRITE (6 , 25  75  ) 

FORMAT ( ///•  • 

STOP 
END 


SET' ,12///// 


IPATH 

25X, '  LU8  EXCEEDED  BY  PATH  #  *,I2) 
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